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Abstract. We have calculated an evolution of protoplanetary disk from an extensive set of initial conditions using a time- 
dependent model capable of simultaneously keeping track of the global evolution of gas and water-ice. A number of sim- 
plifications and idealizations allows for an embodiment of gas-particle coupling, coagulation, sedimentation, and evapora- 
tion/condensation processes. We have shown that, when the evolution of ice is explicitly included, the location of the snowline 
has to be calculated directly as the inner edge of the region where ice is present and not as the radius where disk's temperature 
equals the evaporation temperature of water-ice. The final location of the snowline is set by an interplay between all involved 
processes and is farther from the star than implied by the location of the evaporation temperature radius. The evolution process 
naturally leads to an order of magnitude enhancement in surface density of icy material. 
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1. Introduction 

Water-ice in the protoplanetary disk can exist beyond a mini- 
mum distance from the star called the snowline. Although the 
snowline is defined directly by the presence of ice, the stan- 
dard method of calculating its location is indirect and uses disks 
temperatures. Thus, in the present literature, the snowline is the 
location where the disks temperature is equal to the sublima- 
tion/condensation temperature of water-ice. The total surface 
density of all solids, E s , increases rapidly outside the snowline 
because water-ice, the most abundant species of solid, becomes 
available and its contribution dominates the value of E s . In the 
core accretion - gas capture scenario of giant planet formation 
high values of l, s are necessary to produce solid cores on time 
scales consistent with the presence of a gaseous nebula. Thus, 
at least in such a scenario, the importance of the snowline de- 
rives from a notion that it marks the inner edge of the giant 
planet formation zone. 

In the context of the Solar System formation, the so-called 
"minimum-mass solar nebula" model has been used frequently 
for quantitative analysis. In this steady-state model the temper- 
ature is derived from radiative equilibrium with the solar ra- 
diation field (Hayashi 1981) and the snowline is located at a 
distance R s i = 2.7 AU from the Sun. Just the presence of ice 
exterior to R s i increases solid abundance by a factor of 4 yield- 
ing S s = 6.8 g cirT 2 at the snowline and S, = 2.7 g cirT 2 at 
r — 5 AU (present, and, by assumption, the original location 
of Jupiter). However, even with addition of ice, the values of 
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"E s are still too low for a rapid formation of giant planets cores. 
Lissauer ( 1987 1 calculated that surface density of solids 5-10 
times greater than this given by the Hayashi model are required 
to grow the Jovian core fast enough to accrete the gaseous enve- 
lope before the solar nebula is dispersed. Stevenson & Lunine 
( 1988 ) proposed a mechanism for a further enhancement of the 
abundance of water-ice just exterior to the snowline by diffu- 
sive redistribution of water vapor through the snowline. They 
calculated that such a mechanism can enhance the value of S s 
by 1 to 2 orders of magnitude, making the rapid formation of a 
Jovian core possible at the location of the snowline. However, 
diffusive redistribution cannot explain core formation for other 
giant planets in the Solar System. 

The discovery of ~ 100 giant extrasolar planets (for a re- 
cent review see Bodenheimer & Lin 2002 1 posed formidable 
theoretical problems for the traditional scenario of giant planet 
formation because most exoplanets were found interior to the 
snowline as calculated from the "minimum mass" model, and 
a significant fraction of those (the so-called hot Jupiters) - 
closer than 0. 1 AU, where both high temperature and insuffi- 
cient amount of matter should have prevented their formation. 
It appears that post-formation orbital evolution such as migra- 
tion of a planet under the influence of tidal torque from the disk 
matter (see e.g. Lin et al. 1200 0) or gravitational close encoun- 
ters between planets leading to significant changes of their or- 
bits (see e.g. Rasio & Ford 1996) need to be invoked to explain 
the orbital properties of hot Jupiters. This notwithstanding, a 
revisiting of the snowline concept is also needed to understand 
the original formation of giant exoplanets, and to explain or- 
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bital properties of those exoplanets that seem to escape a sig- 
nificant post-formation orbital evolution. 

Two important developments call for the re-evaluation of 
the snowline concept. First, it becomes clear from astronom- 
ical observations (for a compact review see Beckwith 1994 1 
that protoplanetary disks, of which a solar nebula is thought 
to be a particular example, are active, ever-changing entities 
with limited life-spans that cannot be modeled successfully 
by a steady-state, phenomenological model like the "minimum 
mass" model. Instead, numerical and analytical models based 
on the concept of time-dependent accretion disk theory have 
been developed (Ruden & Pollack WWR Sterzik & Morfill 
119941 Stepinski 1998). In such models different initial condi- 
tions lead to disks of different character and evolutionary his- 
tory, offering the potential to explain observed diversity in exo- 
planets properties. In an evolving disk, the location of the sub- 
limation/condensation radius changes, and it has to be consid- 
ered as a function of time R evap (t). This function is further pa- 
rameterized by initial conditions. 

Second, it has been recognized that the solid particles 
evolve differently to the gas (see e.g. Weidenschilling & Cuzzi 
1 19931 and references therein). The differences between global 
evolution of gaseous and solid components of the protoplane- 
tary disk have been calculated by Stepinski & Valageas ( 1996 
119971 . They have shown that in the evolving protoplanetary 
disk, the sub-disk of water-ice solids decouples over time from 
the gaseous disk leading to significant departures from a con- 
stant solids-to-gas mass ratio The ice sub-disk evolves towards 
a swarm of icy planetesimals characterized by the value of 
Z s considerably higher than at the beginning of the evolution. 
Within such a model the standard method of calculating the lo- 
cation of the snowline is not useful. Because the gas and the 
solids decouple, the snowline location cannot be based on the 
temperature of the gas. Instead, the snowline R s i{t) should be 
calculated directly from its original definition as the inner edge 
of the region where water-ice is present. 

Recently, the location of the snowline has been recalcu- 
lated. Sasselov & Lecar (2000) considered passive and low 
accretion rate models of a protoplanetary disk instead of the 
"minimum mass" model to calculate the location of the snow- 
line using a standard technique of equating R s / with R evap . They 
obtained R s i « 1 AU, significantly closer to a star than a "dis- 
tant" snowline at 2.7 AU as predicted by the "minimum mass" 
model. The focus of the Sasselov & Lecar work was on using a 
model of the disk that is more consistent with observations than 
the "minimum mass" model. Podolak (2003) used a steady- 
state model of an accreting disk to study how the location of 
the snowline depends on accretion rate, ice grain size, and con- 
tamination of ice by other materials. He defines the snowline 
as the location where the rate of ice grain evaporation equals 
the rate of ice condensation from water vapor. This is a more 
accurate, but still an indirect way of calculating the location of 
the snowline. Podolak finds that in the midplane, where most 
grains are to be found, the snowline is nearly independent of 
grain size and composition, but dependent on accretion rate. 

These efforts do not address evolutionary character of pro- 
toplanetary disk, and, in particular, the decoupling of solids 
from the gas. However, the work by Stepinski & Valageas 



( 1997), recently confirmed by more computationally sophisti- 
cated calculations by Weidenschilling (2003) indicates that the 
redistribution of solids is a very important process in the evolu- 
tion of the disk. Due to coagulation and sedimentation the solid 
particles grow on a time scale that is short in comparison to the 
lifetime of the disk. As they grow, they start to drift inward, 
toward the star, and as they grow to ~ 1 km sizes they settle 
into a fixed Keplerian orbit. This ability of the solids to stop 
migrating inwards allows for a seemingly unlikely situation in 
which R s i does not coincide with R evap . The snowline has to 
be calculated directly from its definition by keeping track of 
an evolving distribution of water-ice. Calculating the snowline 
indirectly from gas properties would result in wrongly placing 
the snowline too close to the star. 

In the present communication we demonstrate the differ- 
ence between R s / and R evap for a broad sample of disk models 
analyzed by Kornet et al. ( 120011 f2002i. We identify the radial 
drift as the major factor responsible for the final location of 
the snowline, and we show how redistribution of solids can en- 
hance the abundance of solid material leading naturally to an 
emergence of a giant planet formation zone. The basic assump- 
tions of our disk model and the basic numerical methods em- 
ployed to solve the equations governing its evolution are briefly 
introduced in §2. The results of calculations are presented and 
analyzed in §3, while in §4 we discuss our approach and results 
in relation to ideas and results reported in the literature. 

2. Method of calculation 

The protoplanetary disk is modeled as a two-component fluid 
consisting of gas and solids. The evolution of the gas compo- 
nent is described by an analytic solution to the viscous diffu- 
sion equation, which gives the surface density of the gas as a 
function of radius r and time t (Stepinski 1998). The viscos- 
ity is given by the usual a model. The temperature of the gas 
is calculated in the thin-disk approximation, assuming verti- 
cal thermal balance, according to equations (2) through (6) in 
Stepinski fETO . 

Initially the solids are in the form of grains, but because in 
our model the solids grow all the way to planetesimals our cal- 
culations handle solid objects with a wide range of sizes. We 
refer to them collectively as "dust", although term "solid par- 
ticles" and "dust particles" are also used. The model of their 
evolution includes gas drag effect, sedimentation, coagulation 
and evaporation. Only one component of dust is considered, in 
this paper corresponding to water-ice, which has a sublimation 
temperature T evap — 150 K and a bulk density g - 1 g cm" 3 . The 
main assumptions used in the calculation are (1) at each radius 
the particles are all assumed to have the same size (which, of 
course, varies in time), (2) all collisions between particles lead 
to coagulation, (3) in disk regions with temperature exceed- 
ing Tevap all water is in the form of vapour and evolves at the 
same radial velocity as the gas component, (4) initially, in disk 
regions with temperature below T evap , the particles have sizes 
equal to a mm = 10~ 3 cm, (5) the systematic radial velocity of 
grains is entirely determined by the effects of gas drag, (6) the 
evolution of solids does not affect the evolution of the gaseous 
disk The evolution of solids is governed by the set of two equa- 
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tions (see Kornet et al. (2001 1. The first of them is the standard 
continuity equation for S s . The second equation can be inter- 
preted as the continuity equation for size-weighted surface den- 
sity of solids "L a {r) = a(r)l, s (r) with the source term accounting 
for the growth of solid particles due to collisional aggregation. 
The radius of a solid particle at a distance r from the star is 
denoted by a(r). The set is solved numerically. Details of our 
numerical method is given in Kornet et al. ( 1200 1 i . 

The initial conditions are parameterized by the quantities 
too (the mass of the disk in units of solar mass M ), and jo (the 
total angular momentum of the disk in units of 10 52 g cm 2 s~'). 
Once too and jo are specified, the analytic solution of Stepinski 
(1998) gives the gas surface density E g ,o(/) = £g(M = 0) that 
serves as an intial condition for the evolution of the gaseous 
disk. 

At t = the ratio and the particle radius a are inde- 

pendent of r (values of 0.01 and «„„,, = 10~ 3 cm are adopted, 
respectively). Thus, the initial condition for the evolution of 
solid disk is £5,0(7") = 0.012 g; o(r). We have checked that the 
results are not sensitive to a particular selection of a m i„ by re- 
calculating some models with a mm = 10~ 4 cm. 

3. Results 

Our results are divided into three grids of disk models. Each 
grid groups models with the same constant value of the viscos- 
ity parameter a and contains complete evolutionary informa- 
tion for disks starting from 99 different initial conditions spec- 
ified by 0.02 < too < 0.2 M and 1< j < 25. The grids are 
characterized by a = 10~ 3 , 10~ 2 and 10 _1 , respectively. The 
models are evolved until either (1) the outer edge of the dust 
disk falls within 0. 1 AU, in which case all dust is assumed to 
have accreted onto the star, or (2) the total elapsed time is 10 7 
yr. In the second case, it usually occurs that the surface density 
distribution of solids, £ s (r), converges to a stationary configu- 
ration well before 10 7 yr. 

An example of such convergence is illustrated in Fig. ^ 
where, for a disk with too = 0.02, y'o = 1 and a = 0.1, the sur- 
face density of solids is shown at several selected times as a 
function of the distance r from the star. Initially the outer edge 
of the disk R out and the sublimation limit R evap are located, re- 
spectively, at ~50 AU and ~4 AU. The surface density of solids 
drops from about 0.5 g ctrT 2 at R evap to about 0.1 g cirT 2 at 
R out . Interior to R evap the solids remain in the form of vapour, 
and as such they are tightly coupled to the gas; however at r > 
R evap the solid particles begin to drift with respect to the gas. 
While drifting, they grow at a rate dependent on r (more slowly 
at larger r because of lower densities and collision rates). By t 
= 500 yr enough solid material from the outer disk has drifted 
across the sublimation limit to cause a noticeable depletion of 
solids between r = 4 AU and r = 10 AU. At t = 1000 yr the par- 
ticles arriving from the outer disk to the sublimation limit are 
large enough to fall onto the descending branch of the drift ve- 
locity curve (shown, for example, in Fig. 1 of Weidenschilling 
( 1997 1), and as a result a maximum of S s develops immedi- 
ately exterior to R evap . Due to large viscosity, the gaseous com- 
ponent of the disk evolves so quickly and cools so rapidly that 
at ~ 5 x 10 3 yr R evap begins to decrease, and the inner edge of 
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Fig. 1. The evolution of solids in an a = 0.1 disk with too = 0.02 
and /'o = 1 . Top: surface density of solids at the times indicated 
in the frame (solids that have sublimated but not yet accreted 
onto the star are not taken into account). Bottom: locations of 
sublimation limit (dashed) and snowline (solid) as a function 
of time The snowline is defined as the inner edge of the region 
where ice grains are present. 



the dust disk moves closer to the star. At about the same time 
the solid particles that have arrived at the sublimation limit are 
so large that they practically stop drifting. As their motion is 
essentially Keplerian, one may say that they form a "Keplerian 
barrier" blocking the flow of solids across Revap • Beyond ~ 15 
AU £ s is now markedly lower than in the initial model, as the 
solids originally residing there have moved closer to Revap • An 
outline of the final dust disk is already visible in a region be- 
tween r = 3.5 AU and r — 15 AU. The evolution is completed 
at t = ~ 3 x 10 4 yr, when all solids that have not evaporated are 
collected in that region. E, is now almost uniform, reaching a 
value ~ 4 times larger than at the same locations in the initial 
model. The typical particle radius a is 300 m. At t — ~ 9 x 10 4 
yr a — 1 km - i.e. the solids have evolved into planetesimals. 
At still later times the planetesimals would tend to grow even 
further, but the model is no longer valid because it does not in- 
clude gravitational effects. The snowline (i.e. the inner edge of 
the planetesimals swarm) stabilizes at R s i = 3.5 AU, while the 
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Fig. 2. The decimal logarithm of surface density enhancement £ „ defined by equation Q for models with a = 0.001 (left) and a 
= 0.1 (right). The values of log are plotted on the plane of mo and jo (the initial values of mass and angular momentum of the 
disk). Within the light-gray region at the lower right corner of each frame the solid component has completely accreted onto the 
star. The step-like shape of the boundary of that region results from insufficient resolution of the grid (the locations of our models 
are indicated by dots). 



sublimation limit continues to shrink, and at t ~ 3 x 10 6 yr it is 
equal to -0.5 AU. 

Note that the movement of the sublimation limit is always 
slower than the velocity of the gas, so the vapor does not leak 
to the region beyond R emp . In addition, unlike in the model 
proposed by Stevenson & Lunine (1988 1, in our model there 
is no significant outward diffusion of water vapor through the 
sublimation/condensation radius. This is because the vapor is 
embedded in a dynamically dominant gas and diffuses together 
with it. In viscous accretion disk the net result of diffusion of 
adjunct rings of gas is to produce an inward mass flow at times 
and radii relevant to the sublimation radius. Thus, there is a 
negligible leak of vapor to the region r > R evap . Overall, in 
our model, the disk develops a zone where T < T evap , but ice 
does not exist, as it has been stopped at larger radii by particle 
growth and an associated decay of particle radial drift. 

In the above example the final planetesimals swarm forms 
a well-defined ring ~3.5 AU < r < ~15 AU,with X, in excess 
of 2 g crrT 2 . Whether such a swarm leads to the formation 
of giant planets cores is unclear as "E s is lower than required 
(see Lissauer 1987 1, but the surface density is relatively con- 
stant throughout the entire zone in contrast to the steep decline 
(Lj ~ r~ 3 ^ 2 ) in the "minimum mass" model. Similar structures 
with different sizes and locations are observed in all models in 
which the solids are present at the end of the evolution. In all 
models the ratio of final-to-initial E s is greater than 1 . The latter 
property of the final planetesimal swarms is illustrated in Fig. [2] 
To obtain Fig.|2] the maximum value of the enhancement ratio 

Jl^it = 10 7 ) 

& = — — (1) 

s E s (fl'™ c , t = 0) 



was calculated for each model, and the distribution of £ s was 
plotted on the (mo, jo) plane. In the above formula Y." S MX is the 
maximum density found in a given disk, and R"' ax is the loca- 
tion at which that maximum was found. The distributions are 
shown for a = 10 -3 and a = 10 _1 (the plot for a - 10~ 2 is 
omitted because it does not qualitatively differ from the two 
cases displayed in Fig.[5J. In each frame the gray area indicates 
the region occupied by disks in which all solids are accreted 
onto the star. The step-like shape of its boundary results from 
the poor resolution of the grid (the locations of our models are 
indicated by dots). The enhancement ratio % s increases as jo de- 
creases and niQ increases, i.e. as the initial models become more 
compact and dense. This is because in denser disks the solid 
particles grow to "Keplerian" sizes more quickly, and once the 
Keplerian barrier is formed at R emp the solids that still remain 
at r > Revap are saved from sublimation. At later times the 
distribution of solids can only evolve toward a more compact 
configuration, i.e. % s can only grow. The earlier the Keplerian 
barrier forms, the more solids are locked in the outer disk, the 
further proceeds the "compaction", and the larger may be the 
final value of 

Note that for every mo a critical value jo, c can be cho- 
sen such that disks with jo < jo,c are entirely void of solids. 
The transit from disks with maximum £, to disks in which 
all solids accrete onto the star is very abrupt. This is because 
the Keplerian barrier is not formed at all when the sublimation 
limit falls initially so close to R ou , that there is no time for the 
particles to grow to Keplerian sizes before they arrive at R evap . 

One can also see that in more viscous disks % s is larger. 
This is because the particles in hotter disks gain larger drift 
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velocities. From eq. (6) in Kornet et al. (2001 1 we can estimate 
the maximum inward drift particle velocity 



v: 



dr 



p\d?)Vk 

p~\dp~)\d~r)V k 
d\np\ C] 
d\nr)V k 

where p is gas density , P - gas pressure, and C s - sound veloc- 
ity. As a result, the outer radius of the planetesimal swarm in 
hotter disks is smaller, and f, may reach larger values through 
more "compaction". 

Our results indicate that the most important event in the 
evolution of the disk occurs when the first "Keplerian" particles 
appear at R evap . Let Rf vap be the value of R evap at the moment 
the particles arriving at the sublimation limit are so large that 
their drift time 

= R evap 



(a(R evap )) 



is equal to 10 6 yr. In the above formula vv(fl) is the drift velocity 
of particles with a radius a, and a{r) is the radius of particles at 
a distance r from the star (recall that in our model at each r the 
particles are all assumed to have the same size). Fig. [5] shows 
how well Rf vap correlates with the inner radius of the final plan- 

etesimal swarm, i.e. with the final location of the snowline, R ,. 

' ' si 

To explain the correlation note that as the disk cools, R evap de- 
creases at the rate R evap . 

As long as particles at R evap have velocities Vd(a{R evap )) > 
Revap they move inward across R evap and are destroyed main- 
taining R s i at R evap . However, at later times, particle at R eva p be- 
comes larger, decouple from the gas and settle into Keplerian 
orbits. This leads to VdifliRevap - Rfvap)) K< Revap, and while 
Revap continues to decrease with time, the solids lack inward 
velocity to follow it. The finite snowline forms approxi- 
mately at the location where the inner movement of R s i stalls, 
that is at Rf vap . In disk's outer zone (r > R^) the inward motion 
of solids continues for some additional time leading to further 
"compaction" of solids. 

Because many exoplanets have been found relatively close 
to their parent stars (see §1), it is natural to identify models 
leading to the formation of giant planets as close to the star as 
possible. Note, however, that these are not necessarly the mod- 
els with the smallest value of R*,. For given values of a and 
mo, the models with maximum value of jo are extended and 
cool and would evolve to have smallest possible values of R -\. 
However, due to their low densities, these are not the disks that 
are expected to produce planetary cores, thus their snowlines, 
although located close to the star, are not of great interest. To 
identify disk models forming giant planets at the closest possi- 
ble distances to the star would require integration of our models 
with the models of giant planet formation, a task that is beyond 
the scope of the present paper. Kornet et al. ( 2002 1 have probed 
this issue working with models similar to those investigated 
here, but with high-temperature silicates instead of water-ice 
being the sole species of solids. They have found that, for given 
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Fig. 3. Correlation between the final location of the snowline 
R sl and the distance Rf vap at which "Keplerian particles" ap- 
pear at the sublimation limit. Points with error bars indicate 
the location of our disk models on the {R s i,Rf vap ) plane. The 
length of the error bars is equal to the local grid spacing in r In 
the lower frame R, is smaller because disks with smaller a are 

St 

cooler. 

values of a and mo > m c /", the disks with minimum (but still 
swarm-producing) values of jo are the best candidates to form a 
close-to-the-star giant planet at r - 2 AU. With water-ice as the 
sole species of solids the results are qualitatively the same, ex- 
cept the inner edge for giant planet formations is at r - 4.4 AU 
and 8.3 AU for a disk characterized by a — 10~ 3 and 10 _1 
respectively. 

4. Discussion 

In this paper we have concentrated on the issue of the snowline 
in a protoplanetary disk in which solids, including water-ice, 
transform to ever larger particles through hierarchical coagu- 
lation. In such a disk the solids, initially in the form of dust 
grains, grow all the way to planetesimals acquiring a signif- 
icant inward radial velocities in the process. However, these 
velocities are lost when the size of solid particles approaches 
planetesimal size, preventing their penetration into the inner 
disk regardless of its temperature. This evolution of solids oc- 
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curs on a time scale that is short in comparison with the disk's 
lifetime. As a result solids decouple from the gas and the loca- 
tion of the snowline cannot be calculated from properties of the 
gaseous disk, such as the gas temperature corresponding to ice 
evaporation. Instead, the snowline must be calculated from its 
direct definition as the minimum radius at which ice exists. 

Using an example of a disk evolving from a set of par- 
ticular initial conditions we have demonstrated that although 
the snowline (R s i) coincides with the location of an ice subli- 
mation/condensation temperature (R emp ) during early stages of 
disk evolution, at later times R s i remains fixed at R J , whereas 
R evap decreases as the gaseous disk cools. Thus, at these later 
times, there exists a zone, R evap < r < R^ p free of ice despite 
having gas temperatures T < T evap . 

This counter-intuitive result can be understood as follows. 
In an accretion disk with the overall inward mass flux, a ra- 
dial zone is constantly replenished by material from beyond its 
outer radius. At the beginning of the disk's evolution the "no- 
ice" zone described in the previous paragraph was too hot to 
support ice. When the temperature in this zone dropped below 
the condensation level, the zone no longer consisted of its orig- 
inal vapor-rich material, instead it consists of material carried 
from an outer disk that has been depleted of ice by decoupling 
of solids from the gas. 

Because we are considering an evolving disk, the evapo- 
ration radius is a function of time, R eva p(f), and calculating a 
unique snowline using the traditional method of solving equa- 
tion T(r, t) = T evap is not possible. The so defined "snowline" 
would be a function of time and could be located at arbitrarily 
small radii. However, R s [ calculated directly from ice presence 
converges to a single value R^, and constitutes a reasonable es- 
timate for an inner edge of giant planet formation zone. The 
location of R*. depends on the disk's initial conditions and the 

value of a. We have demonstrated that the value of R[, is set by 
the value of Rf vap , an evaporation radius at the moment when 
solids there are so large that they settle into Keplerian orbits. 
Thus, the snowline is determined by a complicated interplay 
between coagulation, sedimentation and gas properties. 

We have shown that in the potential giant planet formation 
zone, r > R there is a significant enhancement of density of 
solids above its initial value. This enrichment is due to drift 
of solids from an outer disk. All solids are "compacted" into a 
ring between R*. and some outer radius that, however, is signif- 
icantly smaller than the outer radius of the gaseous disk. Thus, 
in our model, the relatively high surface density of ice can be 
achieved naturally, as a result of solids decoupling from the 
gas. 

The above conclusions are based on a simplified descrip- 
tion of the processes governing the evolution of protoplanetary 
disks, and one may wonder how robust they are. Below we crit- 
ically assess several assumptions on which our model is based. 

The most radical is the assumption about the size dis- 
tribution of solid particles, which at each distance from the 
star are assumed all to have the same diameter. This assump- 
tion is responsible for formation of a "totally impenetrable" 
barrier at R J ,. However, detailed work by Morfill ( 1985) and 
Weidenschilling J1997b showed that the size distribution of 



solids in a protoplanetry disk quickly converges to a stage in 
which most of the mass is concentrated in a narrow range of 
sizes approching the maximum size. Thus, our approximation 
can be regarded as a reasonable idealization. The emergence 
of an impenetrable barrier in our models is corroborated by 
the behavior of solids in the two-dimensional disk model by 
Weidenschilling ( 2003 1 where the incoming mass tends to pile 
up at a distance where large bodies form, producing a sharp 
transition in both surface density and mean size. We expect 
that, in reality, the barrier at R sl would leak some ice particles 
inward, but they would have a negligible mass. 

Our second assumption concerns the 100% efficiency of co- 
agulation. In reality, one hardly expects that collisions between 
solid ice particles always lead to sticking without fragmen- 
tation. Interestingly, our calculations suggest that if solids in 
the protoplanetary disk indeed transform themselves into larger 
sizes by hierarchical coagulation, and the efficiency of coagu- 
lation is low, the disk would be depleted of most of its solids, 
diminishing the opportunity for planet formation. 

To illustrate how the evolution of solids depends on the ef- 
ficiency of coagulation we have considered again the partic- 
ular model described in §3, but with an additional parameter 
< e < 1 regulating the efficiency of coagulation. To vary 
this efficiency, we multiply by e the source function for particle 
growth by coagulation, /, (see Kornet et. al.. 12001b . The model 
shown in Fig.^was obtained for e = 1. Decreasing e to 0.5 
and 0.25 caused the final outer radius of the icy planetesimals 
swarm to shrink from 15 AU to 7.5 AU and 2.75 AU, respec- 
tively. The final inner radius moved from 2.65 AU to 2.50 AU 
and 2.33 AU, respectively. For e = 0.2 all solids were accreted 
onto the star. This demonstrates how low coagulation efficiency 
inhibits development of giant planets formation zone, at least 
within a paradigm of hierarchical coagulation. 

The inferred location of icy planetesimals zone in the early 
Solar System is between 3 AU and ~50 AU. The estimation of 
an inner edge (i.e. the snowline) is based on dynamical consid- 
erations; a closer snowline could produce Jupiter at close loca- 
tion, which would prevent formation of terrestrial planets. The 
estimate of the outer radius is based on an observation that the 
Kuiper Belt seems to be truncated at ~ 50 AU (Weidenschilling 
2003 ). Of all of our models, the model with a = 10~ 3 , otq = 0.2 
Mq, jo = 9 and e — 1 produces the icy planetesimal swarm 
closest to the one inferred for the Solar System. The evolution 
of solids surface density in such a model is shown in Fig. 
the final swarm extends from 3 AU to ~ 25 AU with surface 
density ranging from ~ 8 to ~ 12 g ctrT 2 . Note that the solid 
surface density at 5 AU is above 10 g ctrT 2 , which was the stan- 
dard value used by Pollack et al. J1996i to form Jupiter in less 
than 10 Myr. 

Our model provides a reasonable scenario for the formation 
of icy planetesimal swarms in the Solar System. Note, however, 
that a good fit requires e m 1 . For e = 0.5 and e = 0.25 the outer 
radius of the swarm decreases to 15 AU and 10 AU, respec- 
tively, not far enough to account for the Solar System. Indeed, 
to account for the entire Kuiper belt, models with e increas- 
ing with the distance from the star, up to values greater than 
unity in the outer disk are needed. Note that given the numer- 
ous approximations on which the form of / is based, e > 1 is 
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Fig. 4. The evolution of solids in an a = 0.001 disk with mo 
= 0.2 and jo = 9. Top: surface density of solids at the times 
indicated in the frame. Bottom: final distributions of solids ob- 
tained for the indicated vaules of the coagulation efficiency pa- 
rameter e. 

not as incongruous as it appears. For example, in low-viscosity 
disks (a <~ 0.01) the relative particle velocities induced by 
turbulence become smaller than the differential drift velocities 
(Weidenschilling, 1997). The latter are not included in /, and, 
as a result, in such a regime our approach underestimates the 
particle growth rate which could be countered by increasing 
the value of e. Also, in our models the coagulation rate is un- 
derestimated by the absence of other species of solids. Their 
presence would increase the solid density and thus coagulation 
rate. This also can be qualitatively countered by increasing the 
value of e. 

These factors notwithstanding, a high coagulation effi- 
ciency is needed in our model to produce an icy planetesimal 
swarm leading to a Solar-like planetary system. However, the 
high efficiency of coagulation is questionable. Indeed, a so- 
called "meter-sized particle barrier", wherein particles growing 
to meter-size rubble achieve high relative velocities and their 
collisions lead to fragmentation, is discussed in the literature 
(see, for example Weidenschilling & Cuzzi [19931 . Thus, the 
demand for a high coagulation efficiency in our model, to pro- 
duce planets, can be view as an argument against a hierarchi- 



cal coagulation as the sole agent of solid aggregation. Perhaps 
some form of a collective process to transform solids from dust 
directly into planetesimals, thus bypassing the interim stages 
and avoiding large drift velocities that lead to the loss of solids, 
has to take place for the disk to create planets. Such a mecha- 
nism, in the form of gravitational instability of the dust layer, 
was originally proposed by Goldreich & Ward ( I1973> . but re- 
jected when it became clear that the presence of turbulence 
would prevent dust from settling into a thin enough layer to 
become gravitationally unstable. Later, a variant of the gravi- 
tational instability mechanism was proposed wherein the dust 
density was enhanced in gaseous vortices to produce clumps 
dense enough to be collapsed into planetesimals by gravita- 
tional instability (Tanga et al. 1199 6 1. Barring the existence of 
some collective process, the paradigm of hierarchical coagu- 
lation requires a high efficiency of coagulation for the forma- 
tion of giant planets. This requirement is independent of the 
details of our model and is rooted in the high drift velocities of 
intermediate-size particles. 

Finally, we stress that the models discussed in the present 
communication are not meant to be quantitative. Instead, they 
should be viewed as an illustration of the basic processes in- 
volved, which, we believe, is qualitatively correct despite its 
simplicity. Future improvements of the model should allow for 
a range of particle sizes at each position in the disk and for in- 
teractions between different types of particles; a more detailed 
description of the gas component is also a desirable option. 
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